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The modified O'Raifeartaigh model from the context of the generalized seesaw mechanism of 
neutrino mass is investigated. In our evaluation of effective potentials of the theory, both the 
component field and the superspace formalisms to approach the problem are presented. In the 
component field formalism, we take into account the Bose-Einstein condensates in the scalar sector 
by the method of many-boson theory, i.e. we consider both the condensates and the Hartree-Fock- 
Bogoliubov-type self-energies of quantum fluctuations. The diagonalization of the mass matrix 
", of the fermion sector gives the same functional forms of the mass eigenvalues in the generalized 

f*^ ■ seesaw mechanism. The stability condition in the vicinity of the classical vacuum which shows 

' the generalized seesaw situation is obtained by the examination of the mass eigenvalues of the 

. scalar sector of the model. The superspace formalism will be devoted to a comparison between its 

result with that of the component field formalism. (Keywords: Supersymmetric Effective Theories, 
Neutrino Physics, Supersymmetry Breaking, Nonperturbative Effects) 
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The standard model of elementary particles is the most important achievement in modern physics, and still it 
Qh| gives us the horizon of particle phenomenology [1]. Supersymmetry ( SUSY ) [2-22] can be understood as one of the 
ways toward " beyond the standard model" , from a viewpoint of particle phenomenology. In such a supersymmetric 
approach, a theory has both fermionic and bosonic degrees of freedom, and they interact with each other under a 
supersymmetric manner. Dynamics of interacting boson gases quite often show Bose-Einstein condensation ( BEC ) 
as a universal phenomenon. BEC was first found by Bose [23] and Einstein [24] , and theory of BEC of an interacting 
nonrelativistic boson gas was constructed firstly by Bogoliubov [25]. The Bogoliubov theory has a very universal 
character, and it is the case that the theory can be applied to various interacting boson systems. Moreover, methods 
and concepts of BEC and superfluidity of boson gases are useful to examine/understand an interacting fermion system, 
i for example, BCS ( Bardeen-Cooper-Schrieffer ) superconductivity [26] or chiral condensations in the Nambu— Jona- 
Lasinio model ( NJL ) [27,28] and quantum chromodynamics ( QCD ) [29], from the context of spontaneous symmetry 
breakings. On the other hand, a SUSY multiplet must be broken in a model for phenomenology because we have 
not yet found any superpartner. Due to the ( perturbative ) nonrenormalization theorem, a SUSY breaking cannot 
take place in a perturbation theory, and it should be realized in a nonperturbative manner, i.e. a spontaneous SUSY 
L" ' breaking. A lot of modern particle theoreticians consider that a dynamical symmetry breaking is phcnomcnologically 
' prefered for a SUSY breakdown [3,4,5,7,9,10,12]. 

X ■ 

The most important problem in modern particle physics which has been found by experimental results is the origin of 
masses, their hierarchy, and flavor violations of particles. Recent experimental observations confirmed that neutrinos 
should have very tiny masses, and the seesaw mechanism is one of candidates for providing an explanation to neutrino 
masses [30-34]. Hence, it is an interesting issue to make a SUSY model which will show a seesaw machanism. In the 
ordinary seesaw mechanism, neutrino has both a Dirac and a right-handed Majorana mass terms. The references (35) 
and (36) discussed a generalization of the ordinary seesaw mechanism, added a ( very tiny ) left-handed Majorana 
mass, and some interesting results were obtained. It is well-known fact that the O'Raifeartaigh model breaks SUSY 
at its tree level [14]. Recently, the modified O'Raifeartaigh model has been examined in the context of meta-stable 
SUSY breaking [10-13]. Reference (11) gives mass eigenvalues of scalars and spinors: In fact the eigenvalues take quite 
similar structure with that of the generalized seesaw mechanism [35,36]. The purpose of this paper is to examine the 
modified O'Raifeataigh model under the context of the generalized seesaw mechanism of neutrinos. 

This paper is organized as follows. In Sec. II, we invesitigate the generalized seesaw mechanism of the so-called mod- 
ified O'Raifeartaigh model in the component field formalism. It is suitable to employ several many-body-theoretical 
techniques in the component field formalism though it becomes more lengthy than the superspace formalism. While, 
we have to use the notion of component fields at some discussions also in the superspace formalism, especially if 
we wish to take into account BEC in the scalar sector of the theory. After introducing the modified O'Raifeartaigh 
model, we shortly discuss its symmetry property and the classical minimum. We consider it might be possible that 
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the solution of the classical minimum shows the generalized seesaw mass relation. Then, we employ the many-body- 
thcoretical technique to take into account BEC in the scalar sector. By these preparation, the one-loop effective 
potential is calculated, and stability around the classical minimum will be investigated. Possibility of SUSY break- 
down around the classical minimum also be examined. In a one-loop effective potential calculation, the loop expansion 
must converge rapidly enough, and thus the vacuum of a theory should have a semiclassical nature, will not obtain a 
radical modification by possible quantum corrections. This must be the case in our calculation, and thus we consider 
the situation where quantum corrections around a classical minimum are small. For comparison/supplement to the 
result of the component field formalism, a calculation of the one-loop effective potential in the superspace formalism 
is given in Sec. III. The summary and conclusion of this work is presented in Sec. IV. 

We will follow the textbook of Wess and Bagger for the spinor algebra, gamma matrices and metric conventions 
throughout this paper [2]. ( For example, the metric is j]^ = diag(— 1, 1, 1, 1). ) 

II. COMPONENT FIELD FORMALISM 

A. The Classical Solution 

Our starting point is the following Lagrangian of the modified O'Raifeartaigh model [10-13] of three chiral matter 
fields: 



Here, X, <I>± are chiral ( +; right, — ; left ) superfields. We regard $± as neutrino superfields, mo and denote a 
Dirac and a left-handed Majorana mass parameters, respectively. If X takes a ( very large ) VEV compared with mo 
and tol, then the theory may show a seesaw- type situation in the mass matrix eigenvalues of its fermion sector. The 
usual ( ordinary ) seesaw situation will be achieved by tol — > 0. The mass dimensions of / and g become as follows: 
/; [mass] 2 , g; [mass] . We consider the following global U(l)v ( gauge ) and U(1)a ( chiral ) transformations: 

U(l) v : $+ -> e iav $+, $_^e-"^$_, U(1) A : $+ -» e iotA $+, -> e* QA $_, a v ,a A eH. (2) 

The Majorana mass term of the mass parameter mi explicitly breaks both of these global symmetries. We can choose 
the charge of X to keep the coupling term |X<I> + <I>_ invariant under these transformations as 

U(l) v : X -» e~ 2iav X, U(1) A : X -» e~ 2ioiA X. (3) 

The terms fX and pX^ also explicitly break these global U(l) symmetries. The Majorana mass term will breaks 
the U{1)r symmetry under the following charge assignment of the superfields X and <I>± [10-13]: 

U(Y)r : 9 — > e~* aR 6, -> e laR d, X -» e 2iaR X, $+ -» $+, $_ -» e 2Mi ?$_, a R e R. (4) 

Another i?-charge assignment is also possible: 

U(1) R : X -> X, $+^e Mj? $+, -» e ioiR <$>-. (5) 

In this case, the i?-symmetry will be restored at the limit / — > 0. Therefore, the term fX and the Majorana mass 
term are incompatible with respect to the U(1)r symmetry. As a result, there is no global U(l) symmetry in our 
theory. The absence of i?-axion in the modified O'Raifeartaigh model is discussed in Rcfs. [11,12] from the context of 
meta-stable SUSY breaking, and it is phenomenologically favorable. The ordinary O'Raifeartaigh model corresponds 
to the case tjil = m\ = 0, it has an i?-symmetry, and SUSY is broken at the tree level [12,14]. In the ordinary 
O'Raifeartaigh model, the classical solution becomes — <p- — with <fix = arbitrary, and SUSY is spontaneously 
broken in the vacuum. The one-loop effective potential of the O'Raifeartaigh model was calculated in Ref. [15]. In 
that calculation, the degeneracy of vacua is lifted by the one-loop correction, and the origin of the potential becomes 
the only ground state. There is a Z2 symmetry under <I>± — > — $± in (1). The tree level part of the Lagrangian 
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will be obtained from the scalar potential by employing the Euler-Lagrange equations of the auxiliary fields of chiral 
multiplcts: 



V tree [<j>±,<Px] = \F x \ 2 + \F+\ 2 + \F- 



In literature, the classical solution of V tree is given as follows 



+ 



??1 



g4>x<t>+ + m D 4>- 



+ 



m D 4>+ + 2m L (j)- 



(6) 
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(7) 



Here we have assumed that mo, tul, g and / are real valued to obtain the classical minimum. When these parameters 
are real, the classical solution shows the spontaneous Z 2 symmetry breakdown. All of the VEVs of the solution will go 
to infinity at the limit g — ► 0, and g = is a singular point for the solution. ( This expression of the classical minimum 
of our model has a similarity with the classical solution of a Ginzburg-Landau-type <p 4 model which describes the 
low-energy property of the Ising ferromagnet [37]. Hence V tree seems to have a relation with the Ising ferromagnet. ) 
ytree van i s hes a t the classical minimum and the AT = 1 SUSY of this model is unbroken at the classical level [11,12]. 
We will see quantum corrections to the classical solution through the following loop expansion calculation, though we 
mainly investigate a possibiliy of the generalized seesaw situation in the vicinity of the classical solution in this work. 
As mentioned above, the ordinary O'Raifeartaigh model gives (0x)=arbitray at its classical solution with broken 
SUSY, and the one-loop correction gives the unique vacuum as the origin of the potential ( the origin of (6) gives 
a finite energy and not supersymctric ). By introducing the left-handed Majorana mass term in (1), (j) x asmcal has 
obtained the explicit expression given in (7) and then we can discuss the strength of the VEV {4>x} which would give 
a right-handed Majorana mass parameter with some dynamics of our theory: If we set to_l = from the beginning 
of our model, we cannot obtain an explicit expression for (4>x) ( at least at the classical level ). This is crucial in 
the context of this work. The ordinary seesaw mechanism ( the case ttll — ) cannot be considered by the ordinary 
O'Raifeartaigh model. In this paper, we examine (i) when the classical solution can give the generalized seesaw 
situation, (ii) how the vicinity of the classical soution in the one-loop potential is stable and robust, (iii) whether the 
vicinity of the classical solution in the one- loop potential breaks SUSY or not. Usually, SUSY is broken if there is an 
i?-symmetry in a theory, while SUSY will be kept if an i?-symmetry is broken. 

After eliminating the auxiliary fields of X and <J>± and perform integrations of Grassmann coordinates, one finds 
the expression of C in terms of component fields as follows: 



-\m D \\\^ + I'M 2 ) - |s| 2 |<M 2 l<M 2 - 4|™l| 2 |<M 2 - ^H<M 4 - \(fgd> 2 + + fg^) |/| 2 

-{g^m D cft x + 2m L "m) D )$ + (l>- - (grrV D (j)x + 2m\m D )(jy_4> + 



At W 



The phases of mass parameters and 4>x are defined as 

\cj> x \e i6x , m D = \m D \e i9 ° 



»x 



m L = tol e 



7X,VD,VL 



e R. 



(9) 



We can absorb only two of these phases 6x, &d and 6l by a redefinition of fields. Hereafter, we set mo = rrr D and 
rriL — Tn) L by a field redefinition while keeping the phase degree of freedom of (f>x without loss of generality. Later, 
we will observe that mass eigenvalues of scalars and spinors are functions of 6x- In principle, if we take into account 
phase degrees of freedom of <f>± and <f>x, the classical solution of V tree becomes ( fi = f, g^ = g, m) D = mo, m^ L = rriL 
are imposed ) 



± y ~~ \/cos(20 + ) ± v/cos2(26l + ) - 1, 

10-1 = - m <fl^ +l \cos(0+ -0_)± Vcos 2 (fl+ - fl_) - ll , 

cos(0 x + 0+ - 0_) ± Vcos 2 (0 x +0+ - 0- 



ox = 



2m L 
m D \(j>. 



- 1 



= \<fr±\e 



i8± 



(10) 
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However, the phase degrees of freedoms of scalars are chosen by the vanishing conditions of square roots in (10) such 
as + = 0, 7r, = 0, 7r, Ox = 0, 7r ( totally 16 solutions, all are degenerate ). The solutions (7) at the classical level 
are special cases of them. On the contrary, later, we show it is important to take into account the phase degrees of 
freedom Ox , by our examination of particle mass eigenvalues at one-loop level. It is quite difficult to take into account 
phase degrees of freedom of scalar fields and mass parameters in a complete manner in our calculation of the one-loop 
level ( and also, possible renormalization to them ), and thus we will use (7) frequently in our discussion. Due to the 
Hermiticity of our Lagrangian, we have obtained the quartic terms |<7| 2 |</>+| 4 /4 and |g| 2 |0 + | 2 |0x| 2 as positive ( more 
precisely, non- negative ) definite in (8). This fact guarantees the convergence of functional integral of variable (f> + 
in Euclidean region. These quartic interactions are hard-core repulsive interactions at \g\ > 0, give a stability of the 
scalar sector. 



B. Bose-Einstein Condensation 



From the examination at the tree-level of our theory, we speculate that a BEC takes place in the scalar sector of 
the effective potential of (1) also in the one-loop level. To take into account the BEC under an appropriate manner, 
the scalar fields will be divided into the condensates and their fluctuation parts: 



5x = (Px 



(11) 



where the superscript c indicates the condensation parts of the fields. We should mention that the classical solution 
and Bose-Einstein condensates of 4>x and (j)± are different in principle [38], 

^dassical ^ ^ ^classical ^ ^ ( 12 ) 

because the latter include quantum corrections [38]. We assume condensates are space-time independent. Under the 
decomposition of (11), one finds 

i^±i 2 - i^±i 2 +<t>i4>± +4>U c ± + m 2 , 

2 h ^(^ + ) 2 + 20^ + + (^) 2 , 0t2^(^) 2 + 20 c + t ^ + (^ t ) 2 , (13) 

so forth. Consequently, for example, the Dirac and Majorana mass terms of the scalar sector give terms linear in <j>±. 
In fact, mo and mi have a role similar to chemical potential of a nonrclativistic boson theory. The terms linear in the 
fluctuating fields ( and tadpole- type diagrams ) will be dropped from our Lagrangian. This "variational" condition 
corresponds to the Euler-Lagrange equations for condensates [37-42]. The quartic interactions of scalars in C become 

M 2 ,, l4 _ ls |2r 



4 -|0+| 4 = ^- ^+| 4 + 4|0 + | 2 |^| 2 + |^| 4 + 2\4> + \ 2 {4>+^ + $ + 4>%) + 4>+4>+<l> c M + ~M+^ C + 
| 5 | 2 |^| 2 |^+| 2 = |.g| 2 [w x \ 2 W + \ 2 + \4> c x\ 2 \4>+\ 2 + l^l 2 l^| 2 + \4>x\ 2 \4>+? + J> x <t> c JUx + <t> c Ux4> x 4>+ 

+^4> + ^ x + 4> c +4> c x4>U+ + ^ c l\4>+\ 2 h + Px\4>+\ 2 4>x + ^\4>x\ 2 4>+ + <f> c +\4>x\ 2 4>\] ■ (14) 

Here, we have dropped the terms linear in <fi + or cj)x and their Hermitian conjugates. We will employ the Hartree- 
Fock-Bogoliubov ( HFB ) approximation to self-energies coming from the quartic and cubic interactions between 
fluctuations with introducing the following vacuum expectation values: 

Ji(ar) = (4>l(x)4>+(x)), J 2 (x) = (4>+(x)4>+(x)), 4(x) = $ + {x)ft + (x)), K^x) = (4> x (x)4>x(x)), 

K 2 (x) = ($x(x)$+{x)), Kl(x) = (4>l(x)^ x (x)), K 3 (x) = $ x (x)$+(x)), K\{x) = (fi + (x)4> x (x)). (15) 

Here, J\ = (<P + <f>+) is normal, while J2 = (</>+(/>+) and j\ = {$ + $ + ) are anomalous self-energies, similar notions to 
the case of the BCS-Nambu-Gor'kov theory of superconductivity [26,43,44]. The anomalous self-energies indicate a 
breakdown of particle-number non-conservation in the scalar sector, and this is of course an independent phenomenon 
with the particle-number-non-conservation caused by the Majorana mass term of the fermion sector. In nonrelativistic 
theory of BEC, anomalous self-energies are negative quantities. Therefore, one obtaines 

\^+\ 2 h -> 2M+ + Ml, |^+| 2 ^^2J 1 ^ + 4j> + , |0 + | 4 ^4J 1 0^+ + JU+h + M ] +4>1- (16) 
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By the HFB approximation, the cubic interactions of fluctuations will also be dropped from our Lagrangian. Hence 
we get 



^I0 + | 4 - ^ 



(4Ji + 4|0;| 2 )^ + + (4 + {^) 2 )4> + 4>+ + (J 2 + {4>+) 2 )4>l4>l + W + 



2\It It _i_ I J,c |4 



(17) 



and 



l5| 2 |^| 2 |0+| 2 



Wxvm 1 + (\<t>x\ 2 + m<t> + \ 2 + m\ 2 + JiMx\ 



4*% 



K 3 )^ 



'x 



^X 



</>+ + + Kl)^+cj> x + {4> c +( j> x + K 2 )cj> x <j> 



Ml 



(18) 



From the classical solution, we guess |g| 2 |<?!>+| 2 and |g| 2 |03f| 2 take values of ~ 0(|g|°), and thus a matrix elements 
given by a polynomial of them are not small enough to neglect from our Lagrangian. While, we hope < \g\ <C 1 to 
be satisfied for convergence of a perturbative series/diagrams in terms of \g\. This condition could conflict with the 
seesaw condition g(4>x) 3> 3> m^: A realization of the generalized seesaw situation is a non-trivial problem in our 
theory. In the next subsection, we will evaluate the one-loop effective potential of our theory. We hope the potential 
captures the essential feature of quantum dynamics of the system (1) even at the one-loop level. On the other hand, 
we simply drop — f (^xV'+V'+ + 4>+'4'x^P+) + (h.c.) which will give a coupling between the scalar and spinor sectors in 
C. 

Now, we examine the variational condition, namely the vanishing condition of the terms linear in fluctuating scalars. 
From (8), we obtain the linear terms as follows: 



-4>+ \\m D \ 2 (t>^ + pgcj) c + + gm X D 4> c x 4> c l + 2m[m D ^ 



+ 9 



2( \<t>%\ 2 + J 1 +2\4> x \ 2 + 2K 1 f , 4 + 2Kt 



+ 



2 " + 2 

(|m D | 2 + A\m L \ 2 )4> c l + {g1m D (f# + 2m L m{ J )^ t 



f + +K, 



\4>x + kM 



-<t>x 
+h.c. 



grrJ D <i> c + c 



\9\ 2 (Mx 



K. 



(19) 



They must vanish in our treatment of BEC. All terms given above have mass dimension [Mass] 41 . We will employ a 
kind of Popov approximation to our HFB theory [38] , i.e. J 2 = K 2 = K3 = ( all of the anomalous self-energies will 
be dropped ). At the classical solution (7), the vanishing condition of the coefficient function of <p x gives 



Ji 



4 = -— - 



(20) 



This expression would be modified under a renormalization of bare parameters as J\ = — 2f^ ren " 1 / g^ ren \ We confirm 
that the coefficient of vanishes identically at the classical solution (7). From the vanishing condition of the 
coefficient of <f> + at (7), one finds 



Ki = ^ 1 t =. 9 - 2 ( 



2 (rrr D + 2rriLm 3 D 



<?/)• 



(21) 



This expression of K\ takes a large value at the seesaw condition tud S> tul- 

Next, we will introduce several fields of the following definitions for the convenience of our discussion: 

*=(* x ,*m) T , *x = (V>x,^x) T , ^ m = {^mr^ml) T , *Affl = {ip+A+f, f M i = (f,i) T , 
* = (*x,*m), *x = {-ipx,-i>x), ^m = (^mr,^ml), ^mr = {-^+,->P+), ^ml = {-^-,-^-), 

n = (n x ,n M ) T , u x = (4> x , 4> X ) T , u M = (u MR ,u ML f, u MR = (4> + ,4>lf, u ML = {4>-,4>l) T . (22) 
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Here, ipx is a Majorana, ipMR and V'ml are right- and left- handed Majorana fields, respectively. T denotes transpo- 
sition operation of a matrix. The Lagrangian density will be rewritten in the following form by these fields: 

£ = -v tree [(f> c ± ,(f> c x } + ^n B n + ^ri F i&. (23) 



The matrices f2 B and fi F are defined as follows: 

) B = ( ^XX n XM\ a?. .. = ( Q ±+ Q ±~ \ C,B._=(C,B B \ r,B___ = f^ix 



n« — '"XX "XM —I QD ++ qb — Q B <D B \ Q B — ~k x 

" = I ci B o B I ' "mm = I n« OB I ' ll XM — \ ll X+i ll X-)i "MI = I no 

V MX "mm / V "-+ " — / V "-x 



rn$ x n$ M \ n F xx = i<i), n XM = (-^ c + p + - g ^P-,o), 

\ "MX "MM / 



- / -9<f>°+P+ - gH+P- \ qf _ ( i@ - g<P c x P+ - 9^<t>x p - - m D \ (9 A\ 



The definitions 7 s = -y^ 1 -^ 3 and P± = 1=b ^ 7 '' have been used. The Hermiticity = f2 s is satisfied. Entries of 

il B become such that, 



q b = I □ - IslVi + !</>+ 1 2 ) \ ^ f □ - ml - 4m| 



' VA ' _I □-|g| 2 (J 1 + |^| 2 ) 

1 -¥(4 + (0+ t ) 2 )-4 £ □ -m 2 ,-|fflVi+^i + l</>+| 2 + l0xl 2 ) 



'D 



nB = f -\9?{Kl + <^) -<?W C _ - \g\\K 2 + 

x+ " I - \g\ 2 (4 + <t> c Ui) -\g\ 2 {K, + 4> x 4>%) 

qb = ( -g^m D 4>°l \ o B = ( - m D{gHx + 2m L) \ , 9 ,x 

-gm D p + J> + ~ ~ { -m D (g<j> x + 2m L ) J ' ^ J 

All of the off-diagonal elements of £! B are coming from particle- number-non-conserving interactions and/or mean-fields 
of C under the HFB approximation. The diagonalizations of fi B and Sl F will give " quasiparticle" excitation energy 
spectra of scalar and spinor fields in terms of the bare parameters/fields. Especially we have an interest on whether 
the spinor ip x becomes massive or not under the one-loop quantum correction. If there is no massless fcrmion, then 
there is no Nambu-Goldstone ( NG ) fermion, and the Nambu-Goldstone theorem implies the absence of spontaneous 
SUSY breaking in our theory [9,12]. 

C. The One-loop Effective Potential 

In this subsection, we will evaluate and examine the one-loop effective potential of our theory. We obtain the 
generating functional of our theory as follows: 



Z = J V<j>xV0 x V0 + V<f>l_V0-V0lv^xV^xV^+V^+V^-ViJ 

■ J d 4 z{-F tree [</4,«^] + iirf! B n + + ( 



x exp 



source 



(26) 



The one-loop contribution to the effective potential is evaluated to be 

vW[4>° ± ,<p x ] = v B w + v F w, 

V B ^ = ilnDctfi* = ^lnDctfi MM + ^lnDct(of x - fif M -i— n«A 

z z z \ "mm 7 

V F « = -^lnDct^ = -^lnDetO MM - ^lnDctfe - Q F M -i— Q F MX ). (27) 
z z z v "mm 7 

The effective action is found to be 

r (compo) = -i\nZ = Jd 4 x(-V tre *[<t>l,<p x }-VW[<Pl,<p x }). (28) 
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If we perform the path integration of only Hm and *&m, the generating functional becomes 



Z = 



T>(j)x'D<j)x'D'ipx'Dipx exp 



d*3 



(i) 

A/ 



"^X^XX 



n 



XM 



"mm 



"MX 



12 XAf F 



"MX 



MAI 



*x} 



(29) 



where, 



(i) 



M 



V, 



B(l) 



M 



V, 



F(l) 



M 



V, 



5(1) 



M 



2 lnDetfi MM- 



V, 



F(l) 



M 



-Dctf2 MM . 



(30) 



To obtain this expression of Z, we can regard II x and $>x as Grassmann-even and Grassmann-odd source fields in 
the Gaussian integrations of II m and respectively. 

Let us examine the matrix Q F . Since fijJx is the inverse of propagator of massless fermion, a perturbative expansion 
in terms of (0^ x ) _1 for handling Trlnfr suffers from infrared divergences, indicates that the perturbative expansion 
is an unsuitable method for our model, and thus it is forbidden. If fermion tpx remains massless at the one-loop level, 
the determinant of tt F has a zero point at p 2 = 0, must be factorized like det Q F (p) = (p 2 ) 2 (p 2 + (mass) 2 ) 4 . However, 
the direct evaluation of detQ F (p v = 0) from (24), i.e. at the vanishing four- momentum, gives 



detn F { Pl , = 0) 



(4| ff | 4 |^| 4 m 2 ^ 



l) 



(31) 



Hence, there is no massless particle in the fermion sector at <jf\_ ^ 0, and this fact indicates the absence of a spontaneous 
SUSY breakdown in our theory. Note that this fact is globally the case, whole of the functional space of det Sl F . The 
one-loop effective potential of the contribution of tp± will be obtained after the diagonalization of fi£f M in the following 
form: 



V^ {1) = -zTrln(fc - E F ) 2 (k + E F ) 2 (k - E F ) 2 (k + E F f 



(32) 



The energy spectra E± become 



E F (k) 



Ml = 




Ml 2 , 



\g\ 2 W 



x 



2m| T 2 J ( 



D 



\g\ 2 W 



X 



+ ml + \g\\<t> x \m L cos6 x )- (33) 



Here, the phase of <f>x appears in Ml given above. It was shown in Ref. [36] that a one-loop potential is not degenerate 
with the phase 6x in a Nambu— Jona-Lasinio-type dynamical model of the generalized seesaw mechanism. At the 
case |gi 2 |</>5f| 2 ^ TO r> ^ m i ( satisfied under mo S> mi in (7) ), these mass spectra show the generalized seesaw 
mechanism [35,36], M F is light while Ml is heavy. At mu S> thl, they become 



Ml 



V2: 



M l _ 



2m 2 D 



2m\ 



?Wx\ 2 - 



(34) 



Hence, in the generalized seesaw mechanism, the light spinor aquires its mass of 0(rnh) while mu and |</>x| have 
quite minor contributions to it. By taking into account the result (31), we obtain the following mass formula of the 
fermion sector in terms of the bare quantities: 



((M F ) 2 (Ml) 2 (Ml) 



Fn2\2 



(4|5| 4 !^ 



(35) 



The right hand side becomes 16/ 4 m| / at the classical solution, can take a small value at / — > or m L — ► 0. m L = 
is the case of ordinary O'Raifcartaigh model, and in that case Ml = takes place, indicates a breakdown of SUSY. 
Hence, the expression of mass of ?/>x-ficld is found to be 



M F X = 



2\g\ 2 \^ + \ 2 m L 



^/mjy + 4\g\ 2 \(j) x \ 2 m 2 L - 4\g\\<f> x \m 2 D m L cos6 x 



(36) 
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Note that can become imaginary when inside the square root of the denominator of (36) takes a negative value. 
By putting the classical solution (7), one obtaines 

MF = — ~m 9 \m L 

77^2(1 -cos!?*) 

Mx becomes very small and will behave as a pseudo-NG fermion when /, \g\, tul <C tod, and it vanishes at mi = 0, 
while (M^) 2 is always a potitive quantity. It is an interesting fact that M£ will diverge under X — > 0, namely, 
not well-defined in the limit. ( In the case of a non-SUSY dynamical model of the generalized seesaw mechanism, 
6x = 7T is chosen as the vacuum state [36]. ) Therefore, a careful examination on Ox is important ( crucial ) in our 
theory As examined in the subsection A, Ox of the classical solutions will take or tt. M x ~ 0(fgmi J /m 2 - > ) if we 
choose Ox = 7r. We assume there is no spontaneous Lorentz symmetry breaking in the fermion sector. As a result, 
the spectrum of ^-particle must take the following Lorentz symmetric form: 



E F X = y/k 2 + Mp. (38) 
Therefore, we get the one-loop contribution of the fermion sector as 

U F « = ~Trln(fc - E F ) 2 (k Q + E F f(k E F ) 2 (k + E F ) 2 (k a E F ) 2 (k + E F ) 2 . (39) 

We summarize the result of our analysis of the fermion sector: (i) Our model at = corresponds to the ordinary 
O'Raifcartaigh model, has i?-symmetry and SUSY is broken at the ground state, and we have confirmed that an NG 
fermion appears [9,12]. (ii) Our model at uil 0, namely the modified O'Raifeartaigh model [11,12], will give the 
generalized seesaw mechanism at the point (7), while tpx has a finite mass and SUSY seems not broken. We will 
proceed our examination to the boson sector of our theory. 

In our treatment for the scalar sector, first we solve the secular equation det ^f /M = 0. Though the secular equation 
is quartic in the d'Alcmbcrtian □, fortunately, we can diagonalize Oj^- M analytically because its secular equation will 
be factorized into a product of two quadratic equations of □. The results is 

detftlU = (fc -^fn+)(fco+< 1+ )(fco-<i-)(fco+< 1 _) 

x(fc - £f/ 2+ )(fco + E* 2+ )(k - E* 2 _)(k + Ef n _). (40) 

There is no degeneracy in the spectra obtained from det 0^- M = 0. Here, the energy eigenvalues become 



M2±) 2 -< 



*fB /c 2 + C 3 + |C 4 | ~^~T, I ,\2 , ,7 

M M2± = \l o T^\/{c2-c 3 -\c 4 \) +4|ci| 2 



2 

ci = -m D (g ] 4> c ^ + 2m L ), c 2 = m 2 D + 4m|, 



<■:, e mi, + | 5 |Vi+^i + l0+| 2 + l^| 2 ), ^-^(J 2 + (^) 2 )-^- (41) 



| ci 1 2 includes the phase Ox- If we put the expressions at the classical solution for <f> + , J\ and K\ to C3 and C4 with 
employing the Popov approximation J 2 = 0, we get 



M* 1± = M* 2± 



\q\ 2 \<t>y\ 2 „ 9 m'n rn 3 n 5. 



T2 { -ml + ^ + ^^L- -|,|/) + ml (J£UJxL + m 2 + \g\ \f x \m L cos X ) } 



(42) 



Therefore, we conclude 

M F < M* 1± ,M* 2± , (43) 
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at / < 0, g > 0. The mass eigenvalues of M B n+ and M B n+ become tachyonic at 

M^ 1+ ; c 2 (c 3 - |c 4 |) < |ci| 2 , M^ 2+ ; c 2 (c 3 + |c 4 |) < |ci| 2 , (44) 

and an appearance of tachyon indicates the instability of vacuum state [12,45]. In our condition (44) of tachyonic 
masses, c 3 and c 4 include the HFB self-energies. At the classical solution (7) with J 2 = 0, the tachyon condition (44) 
becomes such that 

_5/-i(^ + ^L)>a (45) 

Here we have assumed /, g as real. ( Again, we wish to rewrite that /, Ji, K\, and J 2 have mass dimension [Mass] 2 , 
while g is dimcnsionlcss. ) Hence, the vicinity of the classical solution is stable if / < and g > 0. Due to the HFB 
self-energies and 4> c + , there are several differences between M± and the mass eigenvalues obtained from fl^M'- The 
mass spectra of bosons and fermions are not symmetric ( namely, not supersymmctric ) in our theory. We regard the 
vacuum energy as the order parameter of SUSY-breaking, we must examine the local/global structure of the one-loop 
effective potential to clarify whether the vacuum energy vanishes or not before concluding a breakdown of SUSY. 

In the determinant det{Slxx ~ ^xm(^mm)~ 1 ^mx}> we wish to concentrate upon the vicinity of the classical 
solution (7). At the point (7) with the neglection of J 2 , 0? , in the expression of (25) becomes diagonal. This helps 
us to evaluate the mass eigenvalue of <fix in an analytic manner. Then we get (M B ±) 2 in terms of bare parameters 
as follows: 

(M| ± ) 2 ee \rm(\A\±\B\), 

\(\a\ 2 + (3 2 )b + c\d\ 2 -(ade^+a^e)\ D _ \(0 + fl)(a*b - tfet)| 

A = 1 — m ' B — 1 — m ' 

be — \e\ z be — \e\ z 

a = -|.g| 2 (J 1 + |^| 2 ), 6 ee □ - m 2 , - 4m|, 

c ee □-m 2 ) -|.g| 2 (J 1 +if 1 + |^| 2 + |^| 2 ), d=-gim D <$, e = -m D tf<$ + 2m L ), 

a ee -\g\\Kl + (3 ee -gim D <f> c _ - \g\ 2 (K 2 + cb%^ x ). (46) 

These a,b,c,d,e,a, (3 are matrix elements of fl B ( see (25) ). We also set four- momentum as p u — in the bosonic 
matrices. We examine the stability conditions of M B ±- Especially, we have interest on its behavior in the vicinity of 
the classical solution with the seesaw condition |g| 2 |</>5c | 2 » m 2 D ^> m 2 L . A direct evaluation from (46) with (7), (20), 
(21), and by employing a Popov approximation J 2 = K 2 = K3 = gives 

(M| + ) 2 = (Affj 2 



{m 2 D + 4m 2 L )(3gf+^ + ^) 

48.g 2 / 2 m| + 2g 2 fm 2 D + &gfm 3 D m L + 6gfm% + 3gf^ + • (47) 

mh 2 m L 1 

Hence, if we take into account the seesaw condition mo S> thl, then we obtain g > and / < is the stability 
condition of (M B ±) 2 - A rough estimation gives 

(M| ± ) 2 ~ -2gf~0(gf). (48) 

We conclude that, with taking into account (45) and (48), the generalized seesaw mechanism can take place under 
/ < and g > ( we have assumed that \gf\{mL/m 2 D ) « 1 ). It is worth noticing that (M B ± ) 2 <~ — m 4 /8m| ( 
negative ) if we set K\ = 0, and thus the HFB self-energy K\ is important for stability of the potential at the classical 
solution. Finally one finds 

E B X ± = ^k 2 + (MZ ± ) 2 , (49) 
and we obtain the one-loop contribution of the scalar sector as follows: 

V B ^ = l -Tr(k -E B + )(ko + E x+ )(k -E x _)(k + E B _) 
x(fc - E* 1+ )(k + Eg 1+ )(k - E B ^)(k + 

x(fc - Ef 42+ )(k + < 2+ )(fc - Ef I2 _)(k + E* 2 _). (50) 
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It is a well-known fact that the naive dimensional rcgularization, suitable to keep a gauge invariance in a non- 
SUSY gauge theory, will break SUSY through the regularization. To circumvent of this problem is relatively easier in 
non-gauge models, while the problem is severe in SUSY gauge theories, and the method of "dimensional reduction" 
regularization seems more suitable [18]. Since the Lagrangian we consider here is not a gauge model, here we employ 
a simple cutoff scheme for regularizations of integrals. After performing the four-dimensional momentum integration, 
the one-loop contribution to the effective potential will be obtained as follows: 



vw = 



1 



16tt 2 
+A 4 In 



A 2 {m| 2 + M B2_ + M *\ + + M»\_ + M*% + + M*%_ 2(MP + + Ml 2 )} 
' (1 + M| 2 /A 2 )(l + M| 2 /A 2 )(l + Mg\ + /^)(l + Mf f 2 _/A 2 )(l + Mff 2+ /A 2 )(l + M* 2 2 _/A 2 ) 



-Ml* In(l 



A 2 



M| 2 



■M|i ln(l 



(1 + M| 2 /A 2 ) 2 (l + Mf 2 /A 2 ) 2 (l 
A 2 \ / A 2 



MP/A 2 y 



)+2M-ln(l + J^L) + 2M-ln(l + ^) + 2M-ln(l + 



M| 2 



MP 



-M*\+ ln(l + t^-) - M&!_ ln(l + 

IVI M1+ / 



A 2 



M S2 

M M1- 



^m 2+ ln(l + jj^-) -M^ 2 _ ln(l + ^) 



A 2 



M m 

IV1 M2+ 



M B2 

1V1 M2- 



A 2 \ 
IF 2 / 

(51) 



where, A denotes the four-momentum cutoff. By the standard method for handling the effective potential, namely, 
remove contributions they will vanish at A — > oo ( there is no divergent constant due to Af = 1 SUSY ), one obtains 



1 



16tt 2 



(M| + ) ln(-p) +(M|_) ln(-^J - 2(M$ ln(-^) - 2(m£) ln(-^) - 2(M*) ln(— ) 



+ (M£ 1+ ) 4 m(^)% (M^ln^)^ (Mf f2+ ) 4 ln(^f±)% (Mf f2 _) 4 ln(^f- ) ' 

( M* ± ,M%, M* 1± , M* 2± , Af£ « A ). 

We have arrived at a generalization of the so-called SUSY Coleman- Weinberg potential discussed in Rcfs. [11,12] ( see 
also, Ref. [46] ). In our V™, the one-loop contribution of A-ficld is also included. For obtaining a one-loop potential 
which will not diverge into the negative-energy direction at the limit \(j) x \ — > oo [47], we should impose both 



■Mi 



'Mi 



(52) 



(M£)<(M^ 1± ) 2 ,(M^ 2± ) 



and 



(M|) 2 < (M| ± ) 2 . 
We have known from (43) that (53) is satisfied, while if 

9 X ~ 7T, 1 > g > 0, > f(m L /m 2 D ) > -1, 



(53) 



(54) 



(55) 



then (54) is satisfied. We should set model parameters with respect to these relations. The mass eigenvalues will 
degenerate under using (7) for <\f± with the Popov approximation. Therefore, we will denote them as follows: 



Ml = M* 1+ = M* 2+ , Ml = M* x _ = M* 2 _ 



Since mD,/>mi, especially we have interest on a situation 

Mx <C Mx < M+ < Ml < M F < M B . 
We have arrived at the order of mass eigenvalues and it is the crucial result for our discussion hereafter. 



(56) 



(57) 



Since our effective potential and mass eigenvalues have similarities with those of the Minimal Supersymmetric 
Standard Model ( MSSM ) [19,20], let us utilize some methods/results from it. In theory of MSSM, SUSY is 
cxpricitly broken by a vacuum energy and several soft mass parameters, while it is not broken in the Lagrangian level 
as the starting point of our model. Thus, discussion on renormalization would become simpler than that of MSSM. 
A renormalization-group invariant calculation for renormalization of our 

ytree + y(l) ig gubtle 

because it includes 

many different mass parameters/scales [19,20,48-53]. Since our interest is to examine a possibility of a realization of 
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the generalized seesaw mechanism in the vicinity of the classical solution (7) with taking into account the one-loop 
contribution (52), we concentrate on a VEV of (f> x under the situation (55). Unfortunately, it is difficult to find a 
global minimum of the potential V tree + because it has many parameters, <\f ±1 (j> x , J\ and K\ which should 
be determined variationally. For example, if we put the expressions of classical solution (7) for condensate <p c + to 
reduce variational parameters and try to find a minimum with respect to variation of (j) c x , the potential might give a 
non-vanishing vacuum energy ( hence SUSY is broken ) because this procedure corresponds to a restriction of trial 
functions in the variation: To achieve a true vacuum might be difficult. In the usual prescription of renomalization, 
a running coupling will be used to remove a renormalization point from a theory to get a physical ( renormalization- 
group invariant ) potential, though this procedure is difficult in our case. To make our problem tractable for our 
purpose, we will use the following definition of V^ 1 ) [50,53] by taking into account the Appelquist-Carazzone decoupling 
theorem [54]: 



V(1) = ^EK 2 -(<) 2 )(<) 4ln ^-^ 2 -(Mn 2 )(<) 4 ln^ 2 



/'•'• M 2 

(Z = *,+,-). (58) 

( A quite clear example of the decoupling theorem can be found in Ref. [20]. ) Here, 9{x) is the Heaviside step 
function, it has been introduced to define mass thresholds inside the potential. We have changed the regularization 
method to MS ( modified minimal subtraction scheme ). fj, = e 3 ^ 2 fl, where p, denotes the MS renormalization scale. 
Of course, Mj 8 and M[ are functions of <j) c x . The logarithmic functions appeared in the above equation must satifiy 
|ln(M 2 /^ 2 )| < 1 ( /i: a renomalization point ) for the justification for our loop expansion. We should find the situation 
where U (1) (m) = and ^(V tree + V^) = are simultaneously satisfied [50]. It is a hard task to arrive from the 
complete theory to an effective theory of lowest region in (57) with running parameters with /i. In such a top-down 
approach, as we know from (57), we have totally six decoupling scales until we arrive at the region [i 2 < (M x ) 2 where 
all particles are decoupled, and then V tree alone gives the renormalization-group invariant potential [50,53]. Hence, 
first we wish to consider the case (M x ) 2 < < others. In this case, the potential of the one-loop contribution can 
be written down as follows: 

V (1) = -^(M^ln^. (59) 

Here, we simply have assumed that the effect of decoupled particles is already included by a renormalization of 
parameter. This gives a positive contribution to our one-loop potential. After put the classical solution (7) to <j>± 
of this choose 8 X — n, and take the derivative of V tree + with respect to \4>x\, we get i^x) = m T>/(^9 m L) 
from the stationary condition: We find that the effective field theory of this renormalization point/scale gives the 
same expression for VEV of 4>x with its classical solution, shows the generalized seesaw mechanism. Needless to say, 
we will also obtain {<j> c x ) = m ii/(^9' m L) a t the complete decoupled region /i 2 < (M x ) 2 because = 0. Since the 
improved V tree is the "exact" potential ( with satisfying the matching condition ) [50,53], an observation will find that 
the vacuum is supersymmetric in the energy scale /z 2 < (M x ) 2 . We conclude that, under a reasonable choice of model 
parameters with respect to the conditions for stability of the vicinity of the classical solution (7) and a justfication 
on convergence of the loop expansion, certainly (7) is robust against a quantum correction, will not obtain a radical 
modification, and thus the generalized seesaw mechanism takes place. 



III. SUPERSPACE FORMALISM 



In this section we will calculate the one-loop effective potential in superfield formalism [2,17,55], though we have 
to use components of supcrficlds at several points of our discussions, especially for our consideration on BEC. For 
example, it seems difficult to consider the HFB approximation of quantum fluctuations of scalars in the superfield 
formalism, and thus the self-energies as J\ , K\ , ■ ■ ■ do not appear in our superspace formalism. It is a problem inherent 
in the superspace formalism, and as a result, the one-loop contribution of the superspace formalism is different from 
that of the component field formalism. Moreover, it seems difficult to examine the generalized seesaw mechanism 
by our superfield formalism of one-loop potential because mass eigenvalues of fermion and boson sectors will not 
be derived under a direct manner. Therefore the component field formalism is better to describe the dynamics and 
physical property of the scalar sector with having the BEC. The purpose of this section is to make a comparison 
between the two formalisms. The examination on the generalized seesaw mechanism is beyond scope of this section. 
We employ the background field method, the standard method of superspace formalism [17], to take into account the 
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BEC in our model: 



x = x c + x, x c = ^> c x + eeF x , x = 4> x + e^ x + eeF x , 



(60) 



Similar to the case of component field formalism given in the previous section, again we assume the condensates <jf x 
and <j>± are independent on spacetime coordinates. The Lagrangian will be converted into the following form: 



C = C c + C, 



+ 



(fX c + ?-X c <i> c + <S>l + m D $ c + <f> c _ + m L $L$l) + (h.c.) 



(61) 



2 

Here, we have used the equivalent relations S 2 (6) = —D 2 /4D and 8 2 (0) = —D /4D under the integration of d A x inside 
the action functional of the theory, and have dropped terms linear in the fluctuating superfields X, $±. We have 
introduced several matrix notations defined as follows: 



s = ($+,$_, & + ,&_) T 

M = | 



16 -2 



a ° Jk C * )S 8 (z-z% C^m D ®a l +m L ® l -^ + 9 -X c ® l -±^. (62) 
-°J2-®a J 2 2 2 



The sigma matrices <j v ( the definition: er° = — l2x2> while ct 1 ,^ 2 ,^ 3 are the ordinary Pauli matrices ) act on the 
two-dimensional chirality space (+, — ). The chiral and antichiral delta functions are defined as 



5$±(z') _ D 
1 



\<T,t I l\ 



D 2 

5& ± (z) 4 



The generating functional will be written down in the following form: 



(63) 



= J VXVX^V$>+V& + V$>-V&_ cxp % J d A xC + (sources) 
= J VXVX^cxp J d 4 x C c + (^X^ 



+ -TilnM + g 



9 = ~j d s z J d s z'^J(z)M- 1 6 s (z - z')j\z'). 



(64) 



Here, d 8 z = d 4 xd 2 9d 2 8. To obtain the final expression of Z in (64), we have neglected contributions of (anti)chiral 
sources. The definition of J is 



J 







" g^i^ ® 



By putting the components of X c and we can confirm the fact that C c becomes 

c c = -v tree W ± A c x\- 

Hence the tree level potential is the same in both of the formalisms. Next, we divide M. as follows: 



(65) 



(66) 



2TT 2 



M = M -M', Mo 



D D 



LrD 





2 n! 



® CT 



1 f-^® CT0 





16 



16 

M' = 



1 cr 



■C 



(67) 
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The one-loop effective action |TrlnAl is evaluated to be 

r&per) = ^TrlnM = ^lnDctM + ^Trln(l-Af- 1 M') = ^Trln(l-iM')DAf - 1 



d 8 zln 1- -C f C 



J 2 ^ 
16D 



(68) 



We have dropped | InDetA'fQ 1 because it does not contribute to r^„ per )- _lhc iciauuuo Jvl0 
the commutator [M.q , M.'] = have been used. From the following identity in superspace, 



The relations M^M^ 1 = U~ 1 Mq 1 and 



D 2 D 
16 



5 2 (6» - ^')5 2 (^- 



8=91,9=8' 



the effective potential is found to be 



V, 



(i) 



(super) 



(super) 1 

= 2 



(p 2 + c f c 



1 

2 tr 



A 2 In 



A 2 



J 9 2 9 2 



^l^| 2 ln( 



1 + 



A 



A 2 



\M c f 



A 2 



(A 2 + \M c \ 2 )\M c f 



l^cl 2 = m|,+m| + — rx 



(69) 



(70) 



Of course, the mass dimension of V"/ 1 ^ N is [massl 4 . Both y tree and v/ 1 ^ N vanish simultaneously at the classical 

' [super) L 1 [super) J 

vacuum (7) and SUSY is not broken. To make our calculation on Q in (64) tractable, we approximate M.~ x by 
replacing C — > mp. Then we get 



J d 8 z\g\ 2 \$ c + \ 2 (x - 1 jft + xt ! 



n-m 2 D 



n-m 2 D 



X 



) 



(71) 



Obviously, C/ includes a Kahler potential and (anti)chiral superpotentials of the fluctuating Y-ficld. From a considera- 
tion by the Wick theorem, one finds (□ — m^)" 1 in the Kahler potential corresponds to the propagator {T(j) + (ft + ), while 



mjj /(□ — m 2 ^) in the chiral and antichiral superpotential parts of Q came from (TF + <p + ) and (TFlcjr,), respectively. 
( An examination of mass dimensions of these propagators is also helpful. ) Because (<1>^) 2 A" or (!./(□ — m 2 D ))X are 
chiral superfields, — |q can be inserted between them inside the integration / d 8 z. Therefore we get 



J d 8 zX^X-iQ = d 8 z{X\X)M 



X 



I 2 i 

2 O-r 



4D/n- 



2 l^+l □- 



(72) 



Integration of VXVX^ will give Det 1 Aix, an d this determinant gives a polynomial of F? and F+ . Because 
F1 = F^ = at (7), the one-loop contribution of Det 1 A4x is also vanish and we conclude that SUSY is not broken 
at the classical solution (7). 
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IV. CONCLUSION 



In summary, we have examined the mass spectra of scalars and spinors of the modified O'Raifcartaigh model 
by our evaluation of the one-loop effective potential in the component field formalism, especially in the vicinity of 
the classical solution (7) of the model, from the context of the generalized seesaw mechanism. The BEC in the 
scalar sector has been considered, while the spinor sector has a mathematical similarity with relativistic theory of 
superconductivity [56,57]. Therefore, some parts of our formulation has some similarities with that of theory of 
supersymmetric (color-)superconductivity [58,59], though the intrinsic dynamics of them are quite different. We 
have emphasized that it becomes possible for us to examine the mass spectra for the generalized seesaw mechanism 
of neutrino by introducing the left-handed Majorana mass term ( namely, a modification [11,12] ) to the ordinary 
O'Raifcartaigh model. Our calculation at the one-loop level of the effective potential of the component field formalism 
indicates that SUSY is not broken in the theory due to the absence of an NG fermion, and have confirmed that SUSY 
is not broken at the classical vacuum of the one-loop potential of the superficld formalism. 

In this paper, we have discussed several VEVs of scalars. It is interesting for us to consider some possible relations 
between inflaton of cosmology, scalar fields of (1), a ( generalized ) seesaw mechanism of neutrino, and a spontaneous 
SUSY breaking. The scalar field <fix seems to have a special role in a determination of local/global minima of an 
effective potential of (1), while its VEV determines a right-handed Majorana mass parameter in our theory. Thus 
it is interesting for us to investigate a possible scenario of a relation between <j>x and an inflaton for our further 
investigation. It might be possible to investigate the lightest fermion ipx-fteld could become a candidate of dark 
matter. In the strong CP problem of QCD, a theta angle gives us an important issue on axion. An investigation of a 
relation between 9x and QCD theta angle ( and also the Peccei-Quinn mechanism [60] ) is far beyond scope of this 
paper, though it is also an interesting problem. 
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